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Abstract 

For Bd,s —> ir^K^ and K^*'K^*> decays, the flavor symmetry breaking effects may be particularly small 

pg ■ since the final state interactions should be the same between the corresponding B^ and B s decays due to 

the charge conjugation symmetry of the final states. This is consistent with the newly measured direct CP 

^ \ asymmetry of B s — > -k + K~ . These decays are thus supposed to be important in testing the Standard Model 

00 

CN ' and in probing new physics effects. However, the observation of pure annihilation decay B s — > it + it 



appears to imply a large annihilation scenario with pa ~ 3, in contrast to the case of pa ~ 1 in B u ^ 
decays in the framework of QCD factorization. This seems to indicate unexpectedly large flavor symmetry 
breaking effects between the annihilation amplitudes of B s and B u ^ decays. This apparent contradiction 
could be resolved by noticing that there is a priori no reason to justify the common practice of assuming 



£> ■ the universality of annihilation parameters for different Dirac structures of effective operators. We then 

00 

argue that, for Bj s — > ir^K^ decays, the flavor symmetry breaking effects of annihilation amplitudes have 

^ ' all been included in the initial state decay constants and are thus small. But the flavor symmetry breaking 

effects in Bd s — > K^K^*' decays are likely to be much larger, as part of the annihilation topologies of 
CO 

B s — > KK decay could be related to B s — > tt + tt~ decay. Therefore when new physics effects are searched 

for in these decay channels, care must be taken to consider the potentially large flavor symmetiy breaking 

X' 

H ■ effects in more details. 
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I. INTRODUCTION 

Charmless hadronic B decays, and in particular their CP asymmetries, are very sensitive to new 
physics since the decay amplitudes are either highly Cabibbo-suppressed or loop suppressed in the 
Standard Model (SM). However, it is notoriously difficult to calculate the amplitudes of hadronic B 
decays reliably, due to non-perturbative QCD interactions. These amplitudes are usually evaluated 
using factorization methods, which however are only valid to the leading order of power expansion 
in l/m&. To go beyond the leading power, model dependence may enter. Therefore, in many cases, 
it is hard to distinguish new physics signal from the SM backgrounds. For example, the difference 
of direct CP asymmetries for B° ->■ n~K + and 5+ ->■ n°K + is observed to be -0.126±0. 022 yj], 
which is unexpectedly large since it would vanish in the limit of isospin symmetry. This so-called 
B — > Ktx CP puzzle, as first discovered by the Belle collaboration D2|], might imply new physics in 
the electroweak penguin sector which violates isospin symmetry. However, a mundane explanation 
of large color-suppressed tree amplitude due to non-perturbative QCD is at least equally possible 
(see, for example, D3D and references therein). 

Flavor symmetry is a powerful tool in heavy flavor physics. Generally, SU(3) flavor symmetry 
may receive large corrections at about 20% level, except isospin which is a good symmetry at a 
few percent level. But the flavor symmetry breaking effects could be much smaller in some cases. 
For example, Lipkin [4] noticed that, for B^ s — > ir T K ± decays, the U-spin (d -B- s) symmetry 
breaking effects should be unusually small since the strong phases from final state interactions are 
exactly the same due to the charge conjugation symmetry of the final states. Therefore it could be 
a robust test of the SM vs New Physics to check the relation between the direct CP asymmetries of 
these two decay channels. Interestingly, direct CP asymmetry of B s — > iv + K~ has been measured 
very recently by the LHCb collaboration to be 0.27 ± 0.04 ± 0.01 oD, which is the first observation 
of CP violation in B s decays. This measurement is consistent with the SM relation between direct 
CP asymmetries of Bd )S — > ir T K ± decays. It has also been shown t hat, by a combined use of flavor 
symmetries and factorization method, B ds — > K^'K^*' decays flfjHlJ] may play an important role 
in testing the SM and in probing new physics effects. 

However, as we will see in the following, recent measurements on pure annihilation decays 
B s — > ix + ix~ and B d — ¥ K + K~ may indicate significant violation of flavor symmetry in B ds 
decays, at least for the annihilation amplitudes. The first evidence of B s — > it + it~ decay was 



reported by the CDF collaboration II 1 3M to be 

B(B S -► 7T+7T-) = (0.57 ± 0.15 ± 0.10) x 10" 



(1) 



where the first errors are statistical and the second systematic. It was soon confirmed by the LHCb 
collaboration with 0.37 fb -1 data 111411 as 



B{B S 



IT 7T 



(0.95Zof 7 ±0.13) x 10" 



(2) 



The average of the above measurements gives (0.73±0.14) x 10~ 6 yj. One expects the branching 
ratio of B d — > K + K~ should not be very different from that of B s — > n + n~ as they can be 
related to each other by U-spin symmetry. But the experimental efforts of the CDF and LHCb 
collaborations reveals a surprisingly small result Q1J 



B(B d ->• K + K~) = (0.12 ± 0.06) x 10" 



(3) 



which is several times smaller than the branching ratio of B s — > ix + ix~ and may imply unex- 
pectedly large flavor symmetry breaking effects. As better understanding on the flavor symmetry 
breaking is crucial to separate new physics signal from the SM contributions, we will reinvestigate 
the flavor symmetry breaking effects in annihilation amplitudes, which is important for charmless 
hadronic B decays. 
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The potential importance of weak annihilation amplitudes was noticed first in [|15 - 
charmless B decays and was predicted in perturbative QCD method in [|18l - l20n . Although be- 



2m, weak 



ing formally power suppressed in Aqcd/^6 in QCD factorization method (QCDF) 112 1 - 
annihilation contributions are supposed to be important, together with the chirally-enhanced power 
corrections, to account for the large branching ratios and CP asymmetries of penguin-dominated 
B decays. In soft collinear effective theory [|25h270, it was argued in [|28[ |29|] that annihilation 
contributions are factorizable and may not be significant numerically in charmless B decays. In- 
spired by the experimental progress, there are some theoretical interest 11301 — 13311 recently in these 
pure annihilation decays. The authors of Bill calculated B s — > tt + it~ and B d — > K + K~ decays 
in perturbative QCD method with the results to be in agreement with the experimental data. Refs. 



M 



320 investigated these channels, together with other charmless hadronic B decays, in QCDF 



and found that SU(3) breaking effects should be taken into account for annihilation parameters. In 
addition, a large annihilation scenario seems to be favored in B s decays. While the authors of Q33|1 
discussed the possibility that these decays can also be generated by rescattering from processes 
such as color-favored tree amplitudes. 



FIG. 1. Annihilation contributions at leading order of a s . 

This paper is organized as follows. In the next section, we shall discuss in detail the possible 
flavor dependence of annihilation parameters for charmless B decays to two light pseudoscalar 
mesons in QCDF method. We then conclude with a summary in section III. 

II. ANNIHILATION AMPLITUDES IN QCD FACTORIZATION 

We first briefl y re view the annihilation amplitudes of B — )■ PP decays in QCDF method, one 



may refer to II23L 12411 for the details. The effective Hamiltonian of AB = 1 can be expressed as 



r 10 

H *ff = JEE *«A(aOQ?(aO , (4) 

* i=l p=u,c 

where \ pq = V* b V pq (q = d or s) is a CKM factor, Ci(fj) is the Wilson coefficient which is per- 
turbatively calculable from first principles. The four-quark effective operators Q\ 2 , Ql ...6 an ^ 
Q?,...,io are tree l eve l> QCD penguin and electroweak penguin operators, respectively. These effec- 
tive operators can contribute to the annihilation amplitudes, as shown in Fig. \T\ The basic build 
blocks for pseudoscalar final states may be simplified by taking asymptotic light-cone distribution 
amplitudes and the approximation r£ ~ r£ = r x 

A\ ~ A 2 ~ 27ia s (9(X A - 4 + tt 2 /3) + r 2 x X 2 A ) , A% ~ , 

A{ = A f 2 =0, Al~12Tra s r x {2X 2 A -X A ). (5) 

In the first (last) two diagrams of Fig. [Q the gluon is emitted from the initial (final) quarks. 
Correspondingly, their contributions to the basic building blocks are labeled by the superscript T 
(T). The subscripts 1,2,3 of A\* denote the different Dirac structure of the four-quark operators 
as (V -A)(V-A),(V- A)(V + A) and (S -P)(S+P), respectively. The ratio r x is defined 
by r x = 2m^ / (mb(m q + m s )) with m q the average of the up and down quark masses. X A 
parameterizes the endpoint singularity as 

X A = In mB (l + p A e^A . (6) 

0.5GeVV HA ) 



Notice that <pA is an arbitrary strong phase and normally pa ~ 1 is assumed, which reflects our 
ignorance on the annihilation amplitudes dominated by the soft gluon interaction. In principle, 
Xa may vary not only for different initial and final states but also for different Dirac structure of 
the effective operators. Since different initial and final states can be related to each other by flavor 
symmetry, the annihilation parameters pa and <\>a should only vary mildly for different decay 
channels. However, there is no a priori reason for the annihilation parameters to be the same for 
A\* with different subscript, though in practice they were taken to be universal for simplicity. 
As we shall see in the following, current experimental data may indicate that a universal set of 
annihilation parameters for A\, A % 2 and A% are disfavored. 

It is convenient to define further the annihilation coefficients b's as 



h ° F < 


?iA[ , 


h Cf 


~C 3 A\ + (C 5 + N C C G )A{ 
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h = %C 2 A\ , 
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h cp 

°3,EW ~ J^ 
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h Cf 

OA,EW - -Tp 
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C 10 A\ + C s Al 



(7) 



Numerically, b^^w and b^EW are negligible due to small Wilson coefficients. For the rest of b's, 
only 6 3 contains A^ term, which is actually a dominant one if the magnitude of A\ are not much 
larger than that of A^ . This is a key observation for our later analysis. 



As discussed in (3 



3211 . the recent experimental measurements have revealed that the scenario 
of universal annihilation parameters for all B — > PP decays is in somewhat disagreement with 
pure annihilation decays B s — > n + n~ and B>d — > K + K~ . This point can be seen clearly in Fig. 
[2l where there is no overlap between the regions of annihilation parameters favored by these two 
decays. 



To draw Fig. |2] we have taken the following input parameters 113411 



/b s = 230 MeV , /b = 190 MeV , m b {m b ) = 4.2 GeV , m s {2 GeV) = 95 MeV , (8) 



and the Wolfenstein parameters D35|] 



A = 0.812 



A = 0.2254 



p = 0.144 



fj = 0.342 . 



(9) 



It is then straightforward to evaluate the decay amplitudes of B s — > it + it and B d — > K + K , 
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FIG. 2. Contour plot of the branching ratios of B s — > 7r + 7r~ and B d — > K + K~ decays as functions of the 
annihilation parameters pa and <pA- The solid lines represent the experimental central values and the grey 
regions correspond to one sigma contour. 
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which can be expressed in QCDF as [|24fl 



A(B S -> tt+tT 
A(B d -> K + K~ 



B° V: h V us b 1 + 2h 



U,EW 



B KK [ V ub V ud 



h + 2b 4 + -b±,Ew 






26 4 + -b^Ew 



264 + -bi t EW 



with 



b: 



i—^fBjnfit 



s/2* 



5 i^K = i—pyfBfKJK ■ 



(10) 



(11) 



Considering the SU(3) flavor symmetry breaking, it is not a surprise at all that the scenario 
of universal annihilation parameters does not work. Actually, it has long been assumed (see, for 



example, Refs 



[3 



M 



3811 ) that the annihilation parameters are slightly different between B ud and 



B s decays. Following this assumption, the annihilation parameters have been carefully studied 
in B2I1 . which implied large annihilation corrections of pa ~ 3 for B s decays, in contrast to 
the case of pa — 1 widely used before. Very recently, first observation of direct CP violation in 
charmless Bs decays, A gp{B s — > ir + K~), has been reported by LHCb [|5|] to be 0.27±0.04(s£a£)± 



It has been pointed out in l36ll37ll that a correction factor due to sizable decay width difference in B s system has to 
be included when the experimentally measured branching ratios are compared to the theoretical branching ratios. 
For flavor-specific decay such as B s —> ir + K~, the correction factor is about 1% and can be neglected. But for the 
cases of B s — > tt + it~ , K + K~, the correction factor could be as large as up to 10%. However more experimental 
information such as time-dependent analysis is required to determine the correction factor, so we shall not consider 
this effect in the following. 



B(tf s -»;r'A ) 




FIG. 3. Contour plot of the branching ratios of B s — > -k + tt~, tt + K~, K + K~ decays and direct CP 
asymmetry of B s — > -k + K~ as functions of the annihilation parameters pa and 4>a- The grey (colored) 
regions correspond to one sigma contour. 

O.Ol(syst). The CDF collaboration also reported an evidence of the CP violation to be 0.22 ± 
0.07(stat) ± 0.02(syst) 03911 . A naive average of the latest results yields 0.26 ± 0.04. Including 
this significantly improved data, we confirm that large annihilation scenario for B s — > PP decays 
is still consistent with all the experimental data, as shown in Fig. |3] The overlap regions in the 
figure represent that there exist two solutions satisfying all the constraints: one with pA — 3.5, 
<p A ^ -100° and the other with p A a 3.8, <p A ^ 110°. 

To draw Fig. |3] we have adopted the form factor F B " K = 0.24 as suggested by D38|]. One can 
then easily fit the experimental data of 10 6 B(B S — > ix + K~) = 5.4 ± 0.6q, though it is nontrivial to 
fit the data of 10 6 £(£ fl -> K + K~) = 24.5 ± 1.8. Therefore B(B S -»■ n + K') only gives a rather 
weak constraint on (pa, 4>a) , which is plotted as the light pink region in th e figure . Our results are 



not very sensitive to the parameters of wave functions, so we simply take 11241 . 



40] 



A 



B 



200 MeV 



0.25 



.A 



0.06 



af 



0.25. 



a' 2 = u.Zb , a{- = u.uo , a'i : = U.2t> . (12) 

In a word, the large annihilation scenario for B s decays seems to be good, which is mainly 
required to fit the large B(B S — > n + n~). However, the relatively small B(B d — ¥ K + K~) indicates 
that the case of pa ~ 1 is good enough for B ud decays. This may suggest that the flavor symmetry 
breaking effects are unexpectedly large for annihilation amplitudes between B s and B ud decays. 
2 Unless stated otherwise, we shall always cite Heavy Flavor Averaging Group [l(| for the experimental data. 
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But as pointed out in [|4|], the flavor symmetry breaking effects are probably exceptionally small 
in Bd, s — > n T K ± decays since the final state interactions should be exactly the same for these two 
decay channels. It is then straightforward to obtain a well-known relation between the direct CP 
violations of B ds — > ixK decays 

A CP (B S -> n+K~) B(B Q -> n~K+) r(B s )((m 2 Bs - m 2 K )UF B - K ) 2 

A C p(B° -» 7T-K+) B(B S -> n+K-) r(B°)((m% d - ml)f K F^y ' l } 

which have included explicitly part of U-spin symmetry breaking effects in terms of decay con- 
stants and form factors. Experimentally, left hand side of Eq. (TT3T) equals to 3.0 ± 0.5, while 



42] 



the right hand side equals to 2.2 ± 0.6 if we take the form factor F Bn = 0.26 ± 0.03 Q41 
estimated by light-cone sum rules. So Eq. (fT3T > is consistent within roughly one sigma with the 
experimental measurements and there is no sign of large flavor symmetry breaking beyond decay 
constants and form factors in B ds — > tx t K ± decays. This observation is in a sense disagree with 
the above-discussed scenario of large annihilation magnitude pa > 3 in B s decays together with 
Pa ~ 1 in B Uid decays. This issue is also important in B<i >s — > i^°(*)i^°(*) decays, which may play 
an important role in testing the SM and in probing new physics effects Il6l-ll2ll. It stimulates us to 
reinvestigate the annihilation amplitudes in QCDF and look for different possibility not discussed 
before. 

Notice that, for annihilation part, only b 3 and b^,EW appear in B d>s — > Tt Zf K ± decays. We have 
mentioned before that b 3i EW is negligible due to small electroweak Wilson coefficients and b 3 is in 
general dominated by the A{ term as can be checked from Eq. ©. Since A^ represents the gluon 
emitted from the final quarks, i.e. the last two diagrams of Fig. Q3 the initial state dependence 
of the corresponding annihilation corrections must have all been included in decay constants. 
According to the definition, all the decay constants have been taken outside the building blocks 
A^ . Therefore, A^ is independent of initial state and must be the same for B s and B d decays to the 
same final states. According to the above reasoning, the annihilation amplitudes of B ds — > tv t K ± 
decays are roughly the same and we may take just one set of annihilation parameters (p A , 4>a) 
for these two decays. This is also in agree with the observation that flavor symmetry breaking 
should be very small in B d >a — V n T K ± decays, as noticed first in Q4Q. It is then straightforward 
to determine the annihilation parameters for B ds — > n T K ± decays. The small overlap regions in 
Fig. |4] shows clearly that one set of annihilation parameters are consistent with the experimental 
measurements of the branching ratios and direct CP asymmetries of B ds — > ii T K ± decays. There 
are actually two solutions of (pa, 4>a), one around (1.5, — 50°) and the other around (3.0, 160°). 
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B(tf s -»;r'A ) 




FIG. 4. Contour plot of the branching ratios and direct CP asymmetries of B d s — > tt^K^ decays as 
functions of the annihilation parameters pa and (ft a- The grey (colored) regions correspond to one sigma 
contour with only experimental uncertainties included. 



But concerning that there is no indication yet for large annihilation scenario in B d decays, the 
former solution of (pA, (Pa) = (1-5, — 50°) is preferred. At first glance, this solution seems to 
be in disagreement with the experimental measurements of large pure annihilation decay B(B S — > 
7T + 7r - ). But, a closer look at the annihilation part of the decay amplitudes in QCDF reveals that 
only the basic building blocks A\ 2 appear in B s — > 7r + 7r~ decay, as can be seen from Eqs. (171101 ). 
in contrast to the case of B ds — > K Zf K ± decays which are dominated by the A^ term. It is 
a common practice, just for simplicity, to assume the annihilation parameters pa and (ft a to be 
universal for A\ 2 and A 3 . However, as the subscript of A\ denote different Dirac structures of 
effective operators, there is no a priori reason for the above assumption to be a good approximation. 

Therefore it seems to be reasonable to discuss a scenario with three sets of annihilation pa- 
rameters. Specifically, we introduce the parameters (p A , <ft A ) solely for A^. As A{ is independent 
of the initial state, it is justified to adopt only one set of (p A , 4> A ) for both B ud and B s decays. 
In principle, these annihilation parameters may also vary mildly for different final states, due to 
flavor symmetry breaking effects. But we shall not discuss these effects at this stage. For A\ 2 
terms, however, the experimental data on the branching ratios of B s — > n + n~ and Bd — > K + K~ 
strongly prefer different annihilation parameters between B s and B d decays, as shown in Fig. [2l 
So we shall adopt separately (p l As , <P % a s ) f° r Bs decays and leave (p l A , <ft A ) just for B U:d decays. 



As we have discussed earlier, (p A , (f) A ) can be extracted from B ds — > ti t K ± decays to be 
around (1.5, — 50°). But to fix (p As ,(j) As ), one nas to ^ n ^ besides B(B S — > it + it~), other 
observed B s decays which are also sensitive to the size of A\ 2 . It turns out that currently the 
only candidate is B s — > K + K~ decay. However, its decay amplitude depends not only on A\ 2 
but also on A{. For convenience, we shall simply fix p A = 1.5, <\> A = —50° for A{ to estimate 



the parameters p As and <f) As . The results are shown in Fig. |5(a)| One can see that (p l As , (j) As ) are 
constrained to be either around (2.5, 50°) or around (3.5, —110°). However, the degeneracy can 
not be resolved by the direct CP asymmetry of B s -> K + K~, as A C p(B s -»■ K + K~) ~ -10% 
in both regions. This is because, unlike the branching ratio, Acp(B s — > K + K~) is numerically 
not sensitive to the values of A\ 2 . As for Bd decays, the annihilation parameters (p A , <p A ) are 
only loosely constrained, as shown in Fig. |5(b)| where l0 6 B(B d — > K°K°) = 1.21 ± 0.16 has 
been used. This is mainly because B d — ¥ K + K~ has not been observed yet. Nonetheless, it is 
confirmed that the values of (p A , cf> A ) for B u ^ decays are probably quite different from those of 
{PAs,<f>As) fo * B s decays. 

In short, we confirm in QCDF that the flavor symmetry breaking effects should be exceptionally 
small in B d , s — > 7r =F A" ± decays. This is because the relevant annihilation building block A{, 
as being independent of initial states, should be universal for the same final states. Therefore 
the measurements of these two decays could provide a robust test of the SM vs. New Physics. 
Nevertheless, for B ds — > KK decays, the annihilation building blocks A\ 2 are also involved. 
Unfortunately, A\ 2 could be dependent on the initial state and the experimental data of B s — > 
7r + 7r~ and B d — > K + K~ suggests that the dependence of A\ 2 on the initial state is likely not 
small. That is to say, the flavor symmetry breaking effects may be significant for the annihilation 
amplitudes in B d)S — > KK decays. As B d >s — > K°^K ^ are important channels to probe the 
new physics effects, care must be taken to include the possibly large flavor symmetry breaking 
effects in future studies. 

III. SUMMARY 

A combination of flavor symmetry and QCD factorization has suggested that B d>s — > tx t K ± 
and If (*)lf M decays could be important in testing the Standard Model and in probing new physics 
effects. This is partly because the flavor symmetry breaking effects between the corresponding B d 
and B s decays may be particularly small since the charge conjugation symmetry of the final states 
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(a) (b) 

FIG. 5. Left (right) figure: Contour plot of the branching ratios of B s -)■ K + K~, ir + ir~ (B d -)• K°K°, 
K + K~) decays as functions of the annihilation parameters p l As and cj) As (p A and <fi A ). The grey (colored) 
regions correspond to one sigma contour with only experimental uncertainties included. 

is respected by the final state interactions. Very recently, the first observation of direct CP violation 
in B s decays, Acp(B s — > n + K~), has been reported by the LHCb Collaboration which are well 
consistent with the flavor symmetry expectation in the Standard Model. However, the observation 
of pure annihilation decay B s — > n + ir~ and the surprisingly small results on Bd — ¥ K + K~ 
appears to imply a large annihilation scenario with p A ~ 3 in B s decays, in contrast to the case 
of pa ~ 1 in B u> d decays in the framework of QCD factorization. This seems to indicate large 
flavor symmetry breaking effects between the annihilation amplitudes of B s and B ud decays. In 
QCD factorization, annihilation amplitudes are infrared divergent due to endpoint singularity and 
phenomenological parameterization has to be introduced with model dependence. For simplicity, 
it is a common practice, but without a priori justification, to adopt only one set of annihilation 
parameters for the basic building blocks A\ 2 and A{ corresponding to different Dirac structures 
of effective four-quark operators. We notice that the annihilation amplitudes of B ds — > ir T K ± 
decays are dominated by A^ term, which is independent of the initial state except for the decay 
constant. Therefore the flavor symmetry breaking effects of annihilation amplitudes here must be 
small, in accordance with the corresponding measurements. However, only A\ 2 terms appear in 
pure annihilation decays B s — > n + n~ and Bd — > K + K~. So the experimental results of pure 
annihilation decays strongly indicate a significant dependence of A\ 2 on initial and/or final states. 
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It turns out that all of the building blocks A\ 2 and A^ are involved in B^ s — > K^K^ decays. 
Therefore when new physics effects are searched for in these decay channels, one may need to 
consider the potentially large flavor symmetry breaking effects in more details in the theoretical 
analysis. 
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